Quasi-isotropic cascade in MHD turbulence with mean field 
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We propose a phenomenological theory of incompressible magnetohydrodynamic turbulence in 
the presence of a strong large-scale magnetic field, which establishes a link between the known 
anisotropic models of strong and weak MHD turbulence We argue that the Iroshnikov-Kraichnan 
isotropic cascade develops naturally within the plane perpendicular to the mean field, while oblique- 
parallel cascades with weaker amplitudes can develop, triggered by the perpendicular cascade, with 
a reduced flux resulting from a quasi-resonance condition. The resulting energy spectrum E(ku, k±) 
has the same slope in all directions. The ratio between the extents of the inertial range in the 
parallel and perpendicular directions is equal to 6 rms /Bo- These properties match those found in 
recent 3D MHD simulations with isotropic forcing reported in [R. Grappin and W.-C. Miiller, Phys. 
Rev. E 82, 26406 (2010)]. 
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Introduction - Plasma turbulence in a mean magnetic 
field corresponds to a rather ubiquitous astrophysical set- 
ting. Understanding its nonlinear dynamics which gives 
rise to measurable two-point statistics such as the en- 
ergy spectrum is thus highly important. The current 
state of affairs with regard to the simplified incompress- 
ible magnetohydrodynamic approximation may be sum- 
marized as follows. MHD turbulence with strong mean 
field So mainly forms gradients in directions perpendic- 
ular to the mean field, as the nonlinear cascade is based 
on interactions of Alfven waves propagating in opposite 
directions along the mean field. These interactions are 
resonant, and thus most efficient, perpendicularly to the 
mean field HH3. Predicting the exact form of anisotropy 
and the corresponding scaling of 3D spectra is no simple 
problem. 

The first phenomenology proposed for MHD turbu- 
lence was a weak version of Kolmogorov phenomenol- 
ogy, based on non-resonant coupling between oppositely 
propagating Alfven waves, and ignoring anisotropy [1, [f| 
(IK). It predicted a reduced energy flux and an isotropic 
spectrum with a spectral slope —3/2 instead of the —5/3 
Kolomogorov slope. It has been shown since to hold in 
2D MHD with the rms magnetic field playing the role of 
the mean field 0,0. 

Later on, it was proposed that 3D MHD turbulence 
with guide field would develop along eitherpurely per- 
pendicular directions (weak turbulence, [H Ml) or quasi- 
perpendicular directions (strong turbulence . [loL llll] ) . de- 
pending on the importance of the Alfven wave period 
t a = 1/fcii-Bo compared to the perpendicular nonlinear 
turnover time t^L = 1/kb ~ \/k±b. In the strong case, 
anisotropy is ruled by the so-called critical balance (CB) 
between both times \ = ta/tNL = k±b/k\\Bo — 1, which 



corresponds to quasi-resonant couplings and limits the 
growth of parallel wave numbers during the cascade in 
wave number space, while in the weak turbulence case, 
one has x ^ 1- 

The work presented here is meant to solve two issues. 
First, the two theories of weak turbulence, the isotropic 
one (that applies to 2D, 0, [H) and the anisotropic one 
(3D, 0) are not compatible, as the first one has a spec- 
tral slope —3/2 and the second one a slope —2. It is 
however expected that one should pass from the 3D to 
2D when sufficiently increasing the mean field amplitude 
Bq (see also [HI])- Second, a recent DNS study has re- 
ported a 3D regime with mild anisotropy and uniform 
—3/2 spectral slope in all directions [13| (GM). In this 
regime, the extent of the inertial range remains indepen- 
dent of Bo hi the perpendicular directions, and varies as 
1/Bo in the parallel direction. The purpose of this paper 
is thus to build a phenomenology of 3D MHD turbulence 
which tends to the IK regime in the limit -Bo — > oo, thus 
sharing the properties of the GM regime. 

Some authors 3, H| argued that the GM results are 
an artefact due to measuring the spectrum in a global 
frame, while the true anisotropy can be obtained only by 
measuring spectra in a frame attached to the local mag- 
netic field. While agreeing that the local frame of ref- 
erence is the physical spectral symmetry axis, we think 
that still in the GM regime the properties are identical 
in the local and global frame, basically because the ef- 
fect of frame change amounts to shuffling the symmetry 
axis by an angle of about b rms /Bo = 1/5 which is too 
small to change a quasi-perpendicular spectrum into the 
GM spectrum, and cannot substantially modify the GM 
spectrum itself which fills the whole angular range. 

We propose in this Letter a model of the GM regime 
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combining an IK cascade in the perpendicular plane, cou- 
pled with a slave cascade in the oblique and parallel di- 
rections. In the perpendicular plane, the IK cascade is 
based on the b rms field, which plays the role of the lo- 
cal mean field. Critical balance is assumed to develop 
during a first phase, ensuring a minimum extent of the 
spectrum in the parallel direction. Then, starting with 
the forced parallel and oblique modes not satisfying CB 
(i.e. with relatively large k\\), a slave cascade can be- 
gin. The model describes the evolution of fluctuations 
which do not satisfy the standard CB criterion; it con- 
sists in a flip-flop or ricochet process between the parallel 
and oblique wave vectors directions, each ricochet corre- 
sponding to an increase in either the oblique or parallel 
wave vectors (see below Fig.|3]). A quasi- resonance condi- 
tion, which can be regarded as a generalization of the CB 
criterion, limits the increase at each ricochet. This leads 
to a reduced, IK-like energy flux. An isotropic scaling 
results, and as a corollary, the ratio between the parallel 
and perpendicular inertial range extents is b rms /Bo, as 
obtained numerically in GM. 

Equations and basic models - We start with the 
isotropic Euler case. Formally, the Euler equation can 
be written in Fourier space as 



d t u k = % p ,qku. 



(1) 



It is implicit that in the r.h.s. a sum is made on wave 
vectors p and q such that k = p + q. The equation omits 
to mention vector components and a complicated kernel 
has been replaced by the (dimensionally correct) factor 
k. If only a single triad (k, p, q) is retained in the r.h.s. of 
eq. [TJ no spectrum will form, but taking two nearby tri- 
ads opens the possibility of a cascade. This is the basis of 
the so-called shell models starting with the Desnyansky- 
Novikov equations fl6| in which one uses a discretization 
of wavenumber space into shells with mean radius grow- 
ing in a geometric progression: k n — 2 n ko- A still more 
simplified version of this model amounts to write down a 
flux equation for the energy in each shell of wave num- 
bers, with Ek being the ID energy spectrum: 



u 2 J2 ~ kE k 



(2) 



The variable u n = u(k n ) represents all the modes within 
shell number n. The system reads 

d t ul = F n - F n+1 (3) 

where the energy flux in wave number space is 



F — k ii 6 



(4) 



This system admits a stationary solution with the Kol- 
mogorov scaling u oc fc -1 / 3 and the characteristic cascade 
time ttr(k) of the cascade from k n to fc n +i is 



We now consider the MHD case, in which wave prop- 
agation decreases the coupling along parallel wave num- 
bers, compared to the coupling in the previous isotropic 
case. The MHD equations in Fourier space read for- 
mally, written in terms of the so-called Elsasser variables 



u =F b: 



d t u^ ± ik.B a uf = ^ k=p+q kUpU^ (6) 

This is the analoguous to eq. [TJ Changing variables 

4 = C/ fe ± e ±lfc - Sot (7) 

allows to rewrite eq. |6] in terms of the wave amplitudes 
L/* (using k — p + q) as 

d t U± = S fe=p+9 k U^U^e^- Bot (8) 

(See 0]) for an exact form of the equations in the 2D 
case). In the following we do not distinguish further U + 
from U~ and use u for both fields. 

Quasi-resonance - We start initially with a spectrum 
concentrated at large scales. To examine which triads al- 
low to progress along the Ox mean field axis, we assume 
that only quasi-resonant couplings are non-negligible, 
more precisely that: 



(9) 



where T cor is the correlation time of u p u q , that is T~ l = 
mBJi(pUp,qUq). Assuming that nonlinear coupling are lo- 
cal: 



P ' 



(10) 



and also (which is always true) that the perpendicular 
spectrum is dominant, we find: 

T cor ^ l/(qu q ) (11) 

The quasi-resonance condition (eq. EI) thus reads 



q x B < qu q 
In the following we adopt: 

q x ^ qu q /Bo 



(12) 



(13) 



t^ 1 = F/u 2 = ku 



(5) 



Standard CB phenomenology can be viewed as an appli- 
cation of eq.[13]to all vectors of each triad: k x ~ kuk/Bo. 
A representative triad lies within a plane inclined with re- 
spect to the plane yOz by a small angle ~ Uk/Bo. To be 
able to escape this quasi-perpendicular cascade, we need 
now (i) an initial seed well outside the quasi-resonance 
"cone" q x /q — Uq/Bg (ii) a way to build growing triads 
in the (ku,kx) plane. 

The oblique and parallel cascades We now reason in the 
same way as when dealing with eq. [3] above. There we 
considered an isotropic cascade, assuming a progression 



FIG. 1: The basic triad k = p + q for the parallel (k) and 
oblique (p) cascades, with the quasi-perpendicular wave vec- 
tor (q). 

towards larger wave number in discrete steps from k n to 
k n +i. Here we look for a cascade along the Ox axis but 
with a slower and slower pace, due to the quasi-resonance 
constraint. How much slower? 

To answer, we first consider only perfectly resonant 
triads, hence with q x = (orthogonal triad). Such triads 
do not lead to any progression along Ox. Indeed, the 
triad assumes that the mode u p already exists, with a 
wave vector that has already an x— component equal to 
that of the final wave vector. 

Next we consider quasi-resonant triads as in eq. IT3l 
i.e. with p x — k — q x and q x > 0. Such triads allow 
to progress from k n -\ = k — q x up to k n = k. This is 
illustrated in figs.[T][3] Figs. [l][2] show that the same triad 
intervenes to drive the parallel and the oblique cascades, 
while fig. [3] gives a global view of the cascade by showing 
how the quasi-perpendicular modes allow to progressively 
populate the plane along the parallel and the oblique axis. 

In fig. [3l the angle of the oblique path has been chosen 
to be about tt/4. The important point is that the possi- 
bility to propagate the excitation using the triads shown 
in the figure relies first on the fact that the oblique and 
quasi-perpendicular wave- vectors (reps, p and q) are ini- 
tially excited. Then, as the quasi-perpendicular vector q 
grows, it will necessarily come out of the initial rectangu- 
lar gray region. In fact, the continuation of the oblique 
cascade requires that, as q grows, the parallel component 
q x grows as well. Its maximum value is given by eq. 1131 
For such values of q x to be available, we need that the 
perpendicular spectrum develops (or has developed) ac- 
cording to the critical balance, filling a cone around the 
perpendicular direction, with parallel width varying as 
eq.[H 

Reduced flux in the oblique-parallel cascades - Due to 
the quasi-resonance constraint, only triads with q x /q < 
Uq/ Bq contribute to the cascade along the parallel and 
oblique directions. The main set of triads contributing to 
the cascade in the absence of mean field are characterized 
by < q x < q. Hence we see that the contributing sub- 
set of quasi-resonant triads is a fraction of order q x /q of 



FIG. 2: Sketches of the two basic steps of the oblique cascade 
(left) and parallel cascade (right) from the large scales into 
the unexcited part of the spectral space. Left: From p n to 
p n +i (a) triad in dashed red transfers excitation from p n to 
k n - (b) the (blue solid) triad transfers excitation from k n 
to Pn+i- Right: From k n to k n +i (a) triad in dashed blue 
transfers excitation from k n to p n +i- (b) triad in solid red 
transfers excitation from p n +i to k n +i- 
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FIG. 3: Ricochet process: successive quasi-perpendicular 
wave vectors q„ allowing the coupled cascade between the 
parallel and oblique directions, each step being illustrated 
by a triad illustrated in the previous figure. The quasi- 
perpendicular wave vectors q n and q' n satisfy the quasi- 
resonant condition compatible with critical balance. The 
quasi-perpendicular spectrum is thus supposed to be available 
to drive the coupled parallel-oblique cascades. The shaded 
area represents the initial spectral extent. 

the standard fraction. We thus conclude that the energy 
flux in the oblique/parallel directions is to be reduced, 
compared to the unrestricted isotropic flux with no Bq, 
by a factor R: 

R ~ q x /q ~ u q /B (14) 

Starting from the expression of the energy flux with zero 
mean field written as F — kuj,u q (see eq. 4), we take into 
account the reduction factor R to obtain: 

F = k n u 2 k u 2 q /B (15) 

The expression above for the flux implies that the transfer 
time t\ r {k) for the parallel-oblique spectrum depends on 
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the perpendicular spectrum as tj r (k) 1 = F/u 2 . or 

^(A;)- 1 = ku- x (u-/B ) (16) 

A stationary cascade implies that the flux expression in 
eq. [15] be constant. A consequence of eq. [15] is that the 
sum of the two scaling indices for Uk and u q is —1/2, or 
in other words the sum of the ID slopes of the perpen- 
dicular and parallel-oblique spectra is —3. (See [6j,[l7| for 
the case of z + and z~ spectral slopes in isotropic IK tur- 
bulence). This implies that the scaling will be isotropic 
if and only if the perpendicular slope is 3/2. 

If we remember the scenario developed when com- 
menting Fig. [3] the other basic requirement for the 
oblique/parallel cascades to proceed is that the CB is 
at work, i.e. that the quasi-perpendicular wave numbers 
do have a non-zero parallel component excited satisfying 
the CB. So we have shown that if the CB perpendicular- 
cascade with 3/2 spectral scaling is at work, then oblique 
and parallel cascades with the same scaling are possible. 

The existence of cascades restricted to the perpendic- 
ular direction only with 3/2 spectral scaling is possible 
according to either one of two scenarios, the first one 
being the IK scenario (as in 2D simulations) where non 
resonant-coupling is at work, the second one being re- 
stricted to quasi-resonant coupling but with local scaling 
of the velo city- mag net ic field alignment weakening the 
interactions [111 1 1 81 ] . Re-examining our simulation data 
we found that the local alignment scaling is insufficient to 
explain the 3/2 slope, so it is more likely that the IK sce- 
nario is at work in the GM case, in the 2D perpendicular 
plane. 

The spectral extent, which is controlled by viscosity, 
provides a test of our phenomenology of coupled perpen- 
dicular and oblique cascades. Indeed, at the dissipativc 
scale l/kd, one should find equal transfer and dissipative 
times. This leads to two equalities close to one another, 
depending on whether we consider the parallel or per- 
pendicular cascade: 

uk 2 ~ qu 2 q /B Q ~ ku 2 k /B (17) 

vk 2 ~ qu 2 /b rms ~ kul/b rms (18) 

hence 

kd n /k d± = b rms /B (19) 

This relation has been found numerically in GM. 

We have proposed here a model of the GM regime in 
which the generation of the spectrum in the oblique and 
parallel directions is a genuine effect of the nonlinear cou- 
pling between the oblique and parallel directions. The 
cascade is enslaved to the perpendicular cascade, but 
nevertheless the nonlinear coupling is there. This sce- 
nario might not work for Reduced MHD (RMHD) which 
is often proposed as an equivalent but simpler model for 
MHD with strong guide field. Indeed, in RMHD, any 



increase of the spectrum in the parallel direction is con- 
strained by the CB condition fcii < k±b/Bo. The CB 
condition simply reflects the fact that the parallel gra- 
dients don't develop by their own nonlinear interactions 
but instead are formed by plain linear spatial transport of 
the turbulent signals due to perpendicular eddies. This 
is not so in full MHD where parallel coupling are allowed 
to play their own life, even if limited. 

So, if we start with isotropic conditions with weak tur- 
bulence conditions x C 1, either there is a strict barrier 
to a parallel cascade due to the removal of the nonlin- 
ear parallel coupling as in RMHD, or (in the full MHD 
case), there is no barrier, and the scenario proposed here 
is allowed to develop the GM regiime. 

In conclusion, weak turbulence appears to come in 
three flavors in MHD: (i) the isotropic IK phenomenol- 
ogy which holds in 2D MHD (ii) the strongly anisotropic 
version with no parallel interactions (iii) the mildly 
anisotropic GM regime studied here. The first and third 
versions show fc~ 3 / 2 spectra. The last regime is to be 
found only in MHD. The second regime is probably to be 
found as well in MHD and RMHD. We don't know yet 
what determines the transition between the second and 
third regimes in MHD. This issue is a work in progress. 

We thank G. Belmont, Y. Dong and A. Verdini for 
several fruitful discussions. 
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